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1. Introduction 

During passage through a cluster of galaxies some of the photons of the cosmic mi- 
crowave background (CMB) radiation are scattered by electrons in the hot intracluster 
(IC) gas. The scattering slightly modifies the incident Planck spectrum, imprinting on it 
a unique spectral signature that was first described by Sunyaev & Zeldovich (1972). [The 
impact of hypothetical hot intergalactic gas on the CMB spectrum was calculated earlier 
by Weymann (1966) and Zeldovich & Sunyaev (1969).] This Sunyaev-Zeldovich (S-Z) 
effect is a valuable tool for probing the cluster environment and the global properties 
of the universe. The basic significance of the effect was pointed out in the early work 
of Sunyaev & Zeldovich, and in many papers written during the first decade following 
their original paper. General reviews of the effect and its measurements include those by 
RephaeH (1995a), Birkinshaw (1999), and Carlstrom et al. (2002). 

The major challenge of measuring the effect in nearby clusters was taken up soon after 
the effect was identified, but some decade and a half passed before convincing detections 
were made. Observational results from single-dish radio measurements were reviewed 
by Birkinshaw (1999). Growing realization of the cosmological significance of the effect 
has led to major improvements in observational techniques, and to extensive theoretical 
investigations of its many aspects. The use of interferometric arrays and the substantial 
progress in the development of sensitive radio receivers led to first interferometric images 
of the effect (Jones et al. 1993, Carlstrom et al. 1996). Some 60 clusters have already 
been imaged with the OVRO and BIMA arrays (Carlstrom et al. 2002). 

With the many S-Z projects that will become operational in the very near future, 
measurements of the S-Z effect are about to expand tremendously. The effect is expected 
to be detected in thousands of clusters when planned cluster surveys are conducted. For 
many clusters we expect to have detailed spatial mapping of the effect. Together with the 
much expanded spectral coverage, a large S-Z database will greatly advance our ability 
to fully exploit the potential of using the effect as a precise cosmological probe. The main 
observational challenges in precision S-Z work are discussed in the review of Birkinshaw 
in this volume. As noted, general reviews are available of the S-Z effect and the main 
observational results. Our brief review here is meant to be somewhat more pedagogical. 
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2. - The Effect 

The spectral change resulting from Compton scattering of the CMB by IC gas was 
calculated by Sunyaev & Zeldovich (1972) in the non-relativistic limit based on a solu- 
tion to the Kompaneets equation (which essentially is a diffusion approximation to the 
exact kinetic equation). As was clearly demonstrated by Rephaeli (1995b), the approx- 
imate (though elegantly simple) expressions for the intensity change, A/, obtained by 
Sunyaev & Zeldovich (eqs. 6-9) are not sufficiently accurate at high frc;quencies and 
temperatures. IC gas temperatures span the range 3—15 keV; high electron velocities 
result in relatively large photon energy change in the scattering, requiring a more exact 
relativistic calculation. 

2"1. Total Intensity Change. - Using the exact probability distribution, and a 
relativistically correct form of the electron Maxwellian velocity distribution, Rephaeli 
(1995b) calculated the resulting intensity change in the limit of small optical depth to 
Thomson scattering, r, keeping terms linear in r. In this semi-analytic treatment for the 
calculation of Alt - the change of intensity due to scattering by electrons with a thermal 
velocity distribution - the starting point is the probability of scattering of an incoming 
photon (direction /Uq = cos^o) to the direction /Uq = cos^q is (Chandrasekhar 1950) 



(1) 



/(/^o,/^o) = 



1 



where the subscript refers to the electron rest frame. The resulting frequency shift is 



(2) 



s = In {y' jv) = In 



where i/, v' are the photon frequency before and after the scattering, and (i = v/c ]s 
the dimensionless electron velocity in the CMB frame. It is somewhat more convenient 
to use the variables /z, s and (3 instead of /z, ji' and j3. The probability that a scattering 
results in a frequency shift s is (Wright 1979) 



(3) 



" eV(Mo,Mo) , 



Averaging over a Maxwellian distribution for the electrons yields 



(4) 



V, {s) = 



where 9 ~ kTe/rrieC^. The total change in the photon occupation number along a line 
of sight (los) to the cluster can now be written as 



(5) 



/oo 
[n (xe^) - n (x)] V\ (s) ds 
-oo 



where x is the dimensionless frequency, x = hv/kT, T is the CMB temperature, the 
optical depth is r = ctt / riedl, and ut is the Thomson cross section. The change of 
intensity is simply calculated using AIt{x) = ioX^Ant{x), where io = 2{kT)^ / {hc)"^ . 
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In the nonrelativistic limit the Kompaneets eq. can be readily solved leading to a 
simple analytic expression (Sunyaev & Zeldovich 1972) for A/^, 

(6) Alt = ioyg{x) , 

where the dependence on the cluster gas density (n) and temperature is in the Comp- 
tonization parameter (essentially, an integral over the electron pressure), 

(7) y = j {kTe/mc^)naTdl . 
The spectral function 



(8) 9{x) = 



(e^ - 1)2 



.x(e^ + 1) 
- 1 



is negative for x < 3.83 and positive at larger values of this crossover frequency, 217 
GHz. The magnitude of the relative temperature change due to the thermal effect is 
ATt/T ~ —2y in the R-J region, with y 10~^ along a line of sight (los) through the 
center of a rich cluster. 

The motion of the cluster in the CMB frame induces a kinematic (Doppler) S-Z com- 
ponent, which can be easily calculated in the nonrelativistic limit (Sunyaev & Zeldovich 
1980) 

(9) A/fc = -ioTf3ch{x) , h{x) 



(e^ - 1)2 ' 

where (3c = Vr/c, with Vr the line of sight component of the cluster peculiar velocity. The 
corresponding temperature change for this kinematic component is ATk/T = —t(3c- 

The results of the relativistic and nonrelativistic calculations of AIt{x) are shown in 
fig. 1 for kTe = 5, 10, 15 keV. As is evident, the more accurate relativistic calculation 
yields values that are appreciably different from those based on the nonrelativistic formula. 
The deviations increase with Tg and are particularly large near the crossover frequency, 
where the purely thermal effect vanishes. This frequency shifts to higher values with 
increasing Tg. Also, while the dependence Alt on r is still linear to a high degree of ac- 
curacy, its dependence on Tg is sufficiently non- linear (for typical values of Tg) that the 
intensity change should no longer be determined through the use of the Comptonization 
parameter y. 

While most S-Z observations have traditionally been at (relatively) low frequencies 
(~ 30 GHz, or a; ~ 0.5) where the approximate nonrelativistic description is roughly 
adequate, some 16 clusters were observed also at much higher frequencies (with, e.g. , 
the SuZIE and MITO telescopes, up to x ^ 6.2). In fact, most future S-Z projects will 
observe the effect at several high frequencies. Use of the relativistically exact expressions 
for Alt and Alk is clearly necessary at >> 30 GHz, especially when the effect is used 
to determine precise values of cluster and cosmological parameters. Moreover, since the 
ability to measure peculiar velocities of clusters depends very much on observations very 
close to the crossover frequency, its dependence on Tg, which is approximately given 
(Shimon & RephaeU 2004) in 



(10) 



xo ^ 3.8300(1 + 1.12066* -|- 2.07836'2 - 80.748161^), 
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Fig. 1. - The spectral distribution of AIt/{ioT). The pairs of thick (green) and thin (red) hnes, 
labeled with kT^ = 5, 10, and 15 keV, show the relativistic and nonrelativistic distributions, 
respectively. 



should be taken into account. Note also that because high precision S-Z work en- 
tails use of X-ray derived gas parameters, similarly accurate expressions for the X-ray 
bremsstrahlung emissivity have to be employed (Rephaeli & Yankovitch 1997). In the 
latter paper first order relativistic corrections to the electron velocity distribution and 
electron-electron bremsstrahlung were taken into account in correcting values of the Hub- 
ble constant, Hq, that were previously derived using the nonrelativistic expression for the 
emissivity (see also Hughes & Birkinshaw 1998, and Nozawa et al. 1998b). 

The calculation of Rephaeli (1995b) motivated various generalizations and extensions 
of the relativistic treatment. Challinor & Lasenby (1998) obtained an analytic approx- 
imation to the solution of the relativistically generalized Kompaneets equation, whose 
accuracy was then extended to fifth order in Q (Nozawa et al. 1998a). Sazonov & Sun- 
yaev (1998) and Nozawa et al. (1998b) have extended the relativistic treatment also to 
the kinematic component obtaining (for the first time) the leading cross terms in the 
expression for the combined (thermal and kinematic) intensity change, Alt + A/fe, which 
depend on both and Vr- Since in some rich clusters t ^ 0.02 — 0.03, sufficiently accu- 
rate analytic expansions of the expression for A/ in powers of 8 = kTe/mc^ necessitate 
the (consistent) inclusion also of multiple scatterings of order (Molnar & Birkinshaw 
1999, Itoh et al. 2000, Shimon & Rephaeh 2004). Details on the calculation of the full 
effect are given in these papers. 
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The exact relativistic calculation of A/ does not lead to a simple analytic expression. 
In order to obtain an approximate analytic expression (that can greatly simplify the 
analysis of S-Z measurements), one needs to expand the formal expression for AI in 
powers of (the small quantities) t, &, and /3c- For the resulting expression to be accurate 
to within to ~ 2% for kT^ < 50 keV, Shimon & Rephaeli (2004) included terms up to 
tQ^'^ , t'^Q^, and /3^8'*. [Itoh & Nozawa (2004) have derived a slightly more accurate fit 
to the results of the exact numerical calculation; however, results of this fit are given in 
terms of an array of tabulated coefficients.] Since cluster velocities are expected to be 
generally well below 1000 km/s, terms quadratic in f3c can be ignored. Doing so we can 
write the total intensity change as the sum 



(11) AI/z, = T^/,(a;)e^ + t2 ^/,+8(a;)e^"+i - 



ho{x) + J2fj+i2ix)e' 



where fj{x) = x^Fj{x), with Fj{x) defined in Shimon & Rephaeli (2004). 

2'2. Polarization Components. - Due to the particular angular dependence of 
the cross section, incident unpolarized radiation is linearly polarized if it has a finite 

quadrupole moment. The polarization - which is described in a plane orthogonal to the 
los - is completely specified by the Q and U Stokes parameters. Axial symmetry of the 
scattering (as is evident from the dependence of the frequency change, Equation 2, on 
6, 9' but not on (f) or </)') allows the freedom to select a frame such that U = 0. [Note 
that this is no longer the case in the presence of a magnetic field {e.g. Ohno et al. 2003) 
or when the cluster substructure is spatially resolved (Shimon et al. 2005).] Using the 
above notation Q can be written as 

(12) Q(m) = ^t(1 - Mo) / ^ PMlilJ^'M. 

where P2(/^o) the second Legendre polynomial. From the orthogonality of the Lcgcndrc 
polynomials it is clear that Q depends only on the quadrupole moment of the incident 
radiation as observed in the electron rest frame. 

CMB polarization resulting from Compton scattering in clusters was first considered 
by Sunyaev & Zeldovich (1980). Using Equation (12) they obtained, to first order in r, 

(13) Qi = ^0i^'^^^'^ 



20"^* (e=^ - 1) 



and to second order 



(14) Q. = ^A. 

Both components, which are given here in (temperature) K, depend only on the tangential 
velocity of the cluster, c/3t, and thus their measurement provides complementary velocity 
information to (the leading term in) the kinematic S-Z effect, which depends on the los 
velocity. 

The proportionality of Qi to /3| follows from its dependence on the second derivative 
of the Planck distribution, while Q2 is proportional to its first derivative. This feature 
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Fig. 2. - Patterns of the polarization components described in the text. In cases (a)-(d) /3- 
profile with /3 = 3/2 is assumed, and the cluster peculiar velocity is as indicated by the arrow 
Vt- (a) The oc P^t component, (b) The oc fJtr^ component, (c) The thermal component 
oc {kTe/meC^)T^ at frequencies x < 3.83, and (d) at a; > 3.83. Elliptical cluster with the ratio 
of the principal axes b/a = 0.8 and the symmetry (longer) axis in the X direction (e), the 
symmetry axis inclined at 45° to the picture plane (f), and with the symmetry axis along the 
line of sight (g). Figure adopted from Sazonov & Sunyaev (1999). 



renders these two effects distinguishable through multi-frequency observations. More- 
over, due to the different dependence on the optical depth, the two effects are expected 
to be distinguishable by their different spatial signatures in spherical clusters; these are 
shown in fig. 2a & 2b), along with some of the other polarization components that were 
calculated and described by Sazonov fc Sunyaev (1999) and mentioned briefly below. 
Corrections to the oc component (Equation 13) duo to the random (thermal) compo- 
nent of electron velocity were calculated by Challinor, Ford & Lasenby (2000), and Itoh 
et al. (2000). These corrections generally amount to ~ 10% reduction in Qi- 

Another polarization ciomponcnt is due to multiple scatterings; the first scattering 
induces the requisite quadrupole moment; the radiation is polarized upon second scat- 
tering. This component is proportional to t^6 (figs. 2c & 2d). 

In addition to the above polarization components due to quadrupole moments as- 
sociated with electron motions (both directed and random), polarization arises also 
from the intrinsic CMB quadrupole of the primary anisotropy. This achromatic com- 
ponent is i7) where \/^ is the quadrupole moment of the CMB power spectrum, 



THE SUNYAEV-ZELDOVICH EFFECT 



7 



\fC2 ~ 1.9 X 10~^). The magnitude of this component is typically ~ 20 nK; although 
small, it could possibly be distinguished from the other components by its random distri- 
bution over the sky. It was proposed that the dependence of this polarization component 
on Ci can be used to reduce cosmic variance (Kamionkowski & Loeb 1997), and that the 
evolution of C2 with redshift provides a mean to probe dark energy models (Cooray & 
Baumann 2003). 

In the calculations of the above polarization components it is usually assumed that 
the cluster is spherically symmetric. If the cluster is ellipsoidal the optical depth is 
anisotropic and the polarization patterns are as shown in figs. 2c-2g. 

Measurement of polarized S-Z signals has been a major challenge so far. Based on 
the expected high sensitivity and high spatial resolution of upcoming polarization exper- 
iments, the detection of the major signals is projected to be feasible. This has motivated 
more detailed studies of CMB polarization taking into consideration the complex mor- 
phology of evolving clusters that can be followed in hydrodynamical simulations. For 
example, Diego, Mazzotta & Silk (2003) and Lavaux et al. (2004) explored the polar- 
ization induced by bulk motions during merging of subclusters. The fact that both the 
magnitudes and spatial patterns of polarized cluster signals cannot be realistically ex- 
plored in the context of idealized gas density and temperature distribiitions was made 
clear by the work of Shimon et al. (2005). In their paper the leading polarization com- 
ponents (cx r^0, a and cx t^/3) were calculated in a rich cluster simulated by the 
Enzo code. The maximal polarization levels were found to be a few tenths of a /xK, with 
patterns that are considerably more complex than predicted by the idealized calculations 
whose results are shown in fig. 2. 

In addition to polarization in individual clusters the impact of the cluster population 
on the CMB polarization is of considerable interest. This has been studied by, e.g. Hu 
(2000), and is discussed by Cooray, Baumann & Sigurdson in this volume. 

3. — Power Spectrum & Cluster Counts 

It has been realized early on that the scattering of the CMB in clusters induces spa- 
tial anisotropy (Sunyaev 1977) whose basic properties were already determined in the 
context of a simple model for the distribution of clusters and the evolution of IC gas 
(Rephaeli 1981). More than two decades later it is now well established that this is the 
most important secondary anisotropy on arcminute scales. Since the anisotropy arises 
from scattering of the CMB in the evolving population of clusters, its power spectrum 
and cluster number counts can potentially yield important information on the properties 
of IC gas, the cluster mass function, cosmological evolution of clusters and their gaseous 
contents, as well as some of the global cosmological and large scale structure parame- 
ters. Clearly, therefore, the detailed description of this anisotropy necessitates also the 
additional modeling of gas properties across the evolving population of clusters. 

The traditional approach to the calculation of the S-Z anisotropy is based on the Prcss- 
Schechter model for the cluster mass function, n(M, z), the comoving density of clusters 
of mass M at redshift z. Following collapse and virialization, IC gas is presumed to have 
reached hydrostatic equilibrium at the virial temperature, with a density distribution that 
is commonly assumed to have an isothermal /3 profile. The mass function is normalized 
by specifying the mass variance on a scale of 8/i~^ Mpc, as, a parameter which can be 
determined from the primary CMB power spectrum, from large-scale galaxy surveys, or 
from the observed X-ray temperature function when calibrated by a mass-temperature 
relation (which is currently limited to clusters at relatively low redshifts). The cluster- 
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induced anisotropy has been studied at an increasingly greater degree of sophistication 
and detail (and in a range of cosmological and dark matter models) beginning more 
than a decade ago {e.g. , Makino & Suto 1993, Bartlctt & Silk 1994, Colafrancesco et 
al. 1994). For example, in the latter paper the temperature anisotropy was calculated 
in a flat CDM model including gas evolution. The approach adopted in that paper 
was later (Colafrancesco ct al. 1997) extended to other cosmological models, and to the 
calculation of the mass and redshift distributions of the many thousands of clusters that 
are expected to be detected during the planned Planck survey. The anisotropy and its 
power spectra can also be generated directly from hydrodynamical simulations (e.g. , 
da Silva et al. 2000). And, of course, the range of cosmological models was extended 
to include currently favored ACDM models (beginning with the works of Komatsu & 
Kitayama 1999, Molnar & Birkinshaw 2000, and Cooray ct al. 2000). 

The main features of the power spectra of the anisotropy due to the thermal and 
kinematic S-Z components, and the cluster number counts, were determined already 
in the above mentioned papers. Given the nature of the required input parameters 
characterizing the cosmological and cluster models, it is not surprising that the predicted 
power spectra and number counts span a wide range. We briefly review some of the 
recent (since 2001) results from analytical calculations and hydrodynamical simulations 
in ACDM models. 

Two approaches have been adopted in the calculation of the S-Z induced anisotropy. 

In the analytic approach clusters are described by simple models, which characterize their 
morphologies, temperatures, and their gas content and evolution, usually as functions of 
the virial mass and redshift. For each cluster the profile of the Comptonization parameter 
is calculated and transformed into Fourier space. The overall power spectrum can then be 
quantified by convolving the transformed y-parameter with a cosmological mass function 
characterizing the universal cluster population, traditionally taken to be of the Press & 
Schechter (1974) form. The power spectrum can then be calculated numerically using 
the expression 



where r is the comoving radial distance to a cluster of mass M located at redshift 
z, N{M,z) is the mass function, and Q{M,z) is the angular Fourier transform of the 
AT. In the second approach high resolution dynamical and hydrodynamical simulations 
(incorporating both dark matter and baryonic components) yield large cosmological fields 
which enable the identification of individual clusters and their properties relevant to the 
calculation of y, and the S-Z power spectrum. 

Of course, these approaches have their respective advantages and drawbacks. Simple 
analytic modeling in the first approach can be readily implemented in a computer code, 
and the impact of various input data can be easily assessed. The main disadvantage 
of this approach is the explicit need for detailed modeling of each component of the 
calculation, each with its own uncertainties due to either oversimplified physics, or lack 
of quality data, or both. An example for this is the mass-temperature relation, used to 
attribute a temperature to a cluster with virial mass M . A favored scaling is the virial 
relation 



(15) 




(16) 
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where Tq is the gas temperature of a 10^^ Mq cluster located at redshift z = 0, and 
A{flo,z) is the non-hnear density contrast at redshift z, and h is the value of Hq in 
units of 100 km Mpc^^. In addition to the questionable assumption that chistcrs arc 
relaxed and in hydrostatic equihbrium, there is appreciable uncertainty both in the value 
of To and the scaling of the temperature with redshift. The main current limitations 
of cluster hydrodynamical simulations are the inclusion of a restricted range of physical 
processes, and insufficient spatial resolution. 

As consequence of the redshift independence of the S-Z effect the detection of many 
distant clusters by upcoming projects is very feasible. Obviously, the number of clusters 
that can in principle be detected depends on both their universal population, internal 
properties, and the flux detection limit of the experiment. The number of clusters whose 
S-Z flux exceeds AF^, is 



The lower limit of the mass integral corresponds to the limiting flux from a cluster 
with mass M located at redshift z, and B{M, z) is either unity or zero, depending on 
whether the flux measured from the cluster is higher or lower than the flux limit of the 
experiment. Since the detection of clusters through the S-Z effect cannot reveal any 
information concerning their redshifts, the observable quantity is the cumulative number 
clusters at all redshifts. Note that just as for the power spectrum, theoretical predictions 
of cluster number counts span a wide range that reflects the uncertainties in internal 
properties of clusters and their mass function. 

The power spectrum of the S-Z effect in a ACDM model has been calculated by 
numerous authors (e.g. Komatsu & Kitayama 1999, Molnar & Birkinshaw 2000, Cooray 
et al. 2000, Komatsu & Seljak 2002, Bond et al. 2002, Zhang et al. 2002, Springcl et 
al. 2001) both analytically and from results of simulations. The resulting power spectra 
have essentially a universal shape - a steep monotonic ascent to a peak followed by a sharp 
descent, as shown in the (semi-logarithmic) plot in fig. 4. The curves differ, however, 
in the multipoles at which the peaks fall, the magnitude of the peak power, and - to a 
lesser extent - the width of the peaks, as detailed in the rest of this section. 

The power spectrum was calculated by Molnar & Birkinshaw (1999) for an = 0.2, 
SIa = 0.8, h = 0.5 model, using a Press & Schechter mass function normalized by 
cTs = 1.35, a density fluctuation spectrum with index n = —1.4, and cluster mass range 
bounded by Mmin = IO^/j-^Mq and Mmax = 10^^ Mq. Evolving gas fraction 
was adopted, with the density decreasing according to ~ (1 -|- 2;)"^ at earlier times. An 
isothermal gas was assumed, with a temperature of 8.7 keV ascribed to a local cluster 
of mass 10^^ Mq, and a /3 density profile with (3=1. The resulting power spectrum 
peaks at ^ ~ 2000, with a peak magnitude of £{£ + 1)/2tt ■ Ce ~ IQ-^^ K^. Komatsu & 
Kitayama (1999) employed fl^ = 0.5, Hq = 0.3, and h = 0.7, similarly modeled the gas 
profile and mass function, but adopted as = 1, P = 2/3, Mmin = 5 • 10^'^ Mq, and 
Mmax = 5-10^^ Mq. Moreover, the gas content was assumed to be non-evolving, and 
the core radius was determined by assuming that the gas is in a minimum entropy phase, 
which implies an evolving core size. This model was formulated using a free parameter, 
e, with the values —1, 0, 1. The power spectrum peaks at I ~ 2000 — 4000, with maximal 
power in the range 4—7-10^^^ for the three values of e. Note that the spectral function 
g{x) was divided out, so that in order to compare these results with others calculated 
in the Rayleigh- Jeans region, where g{x) = —2, the power spectrum magnitude must 
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be multiplied by a factor of 4. Cooray et al. (2000) calculated the power spectrum for 
an fio = 0.65, J^a = 0.65 and h = 0.65 model, normalized the Press & Schechter mass 
function with erg = 0.9, and determined the gas profile from a solution to the hydrostatic 
equation with an NFW-distributed dark matter profile. The gas was assumed to be non- 
evolving fraction of the total mass, and isothermal with a (local) 10^^ Mq cluster 
assigned to have a temperature of 5.2 keV. The resulting peak power is ^ 10^^^ at 
£ ~ 2000 - 3000. 

Studies of the S-Z power spectrum based on hydrodynamical simulations include those 
of Springel et al. (2001), and Bond et al.(2002), who adopted the favored SIa = 0.7, 
fio = 0.3, model with erg = 0.9 and index n = 1. The resulting power spectra are broad 
and peak at considerably higher multipoles of up to ^ 10000, with peak power that is 
also somewhat higher than calculated analytically. (Similar results were obtained also 
by Zhang et al. 2002.) 

As an example of how different models may lead to different power predictions we show 
in fig. 3 a comparison of power spectra calculated with three different mass functions, 
Press & Schechter (1974), Lee & Shandarin (1999), Sheth & Tormen (1999), hereafter 
PS, LS and ST, respectively. While the PS and ST mass functions seem to yield approx- 




I 



Fig. 3. - Dependence of the S-Z power spectrum on the choice of mass function. 

imately the same results, the LS mass function generates lower magnitude power at all 
muhipoles by - 30 - 50%. 

From these and numerous other published results it is clear that the predicted power 
spectra span a wide range of values for the width of the peak and its magnitude, reflecting 
mostly (but not only) different assumed values cg = 0.9 and gas properties. The impact 
of a particular choice of parameters on the peak power and its typical scale (£) can be 
readily predicted. Well known is the steep dependence of the power on as (fig. 4), due 
to its location in the exponential part of the mass function; various studies have shown 
that even a mild change of ^ 10% may lead to an order of magnitude difference in the 
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magnitude of the power spectrum. On the other hand, the spectrum is much less sensitive 
to other cosmological parameters, such as Om, implying that the uncertainty in ag is of 
major relevance to precise determination of the power spectrum. The impact of different 
choices of IC gas properties can be easily anticipated; for example, an increase of either 
the normalization of the mass-temperature relation or the central electron density by a 
factor /, would lead to an increase of power by a factor since the power spcctriim is 
proportional to y^. A steeper y profile - as is obtained in polytropic gas models - shifts 
power towards higher multipoles since the effective size of the cluster decreases (Sadeh 
& Rephacli 2004). This is also the case if the gas is assumed not to evolve with redshift, 
since this leads to an increase in the population of distant, smaller clusters. 
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Fig. 4. - Power spectrum predictions for various erg values. Solid, long dashed, short dashed, 
dashed-dotted and dotted curves correspond to crs = 1.32, 1.23, 1.18, 1.00, 0.68, respectively. 



Along with calculations of the power spectrum, many studies have been conducted 
of cluster number counts; the detailed predictions for these too show a large variance. 
Colafrancesco et al. (1997) used a ACDM model (rio = 0.2, I^a = 0.8, h = 0.5, n = 1) and 
four different flux limits (10 — 60 mJy) to assess the distribution of cumulative number 
counts as a function of redshift. Their calculations were carried out within the framework 
of both evolving and non-evolving IC gas scenarios; for a flux limit of 30 mJy they found 
cumulative counts of 2 • 10* clusters to a redshift of z 0.1 in the non-evolving gas 
scenario, and ^ 4 • 10^ clusters when IC evolution is taken into account. In both models 
counts drop drastically around z ~ 1, but with the number of clusters at higher redshifts 
still much larger in the non-evolutionary case. Number counts exhibit the same degree 
of sensitivity to as as the power spectrum; a graphical demonstration of this effect is 
shown in fig. 5. 

In fact, while it is relatively easy to run a code with various input data and examine 
the resulting output, what is of practical interest is to follow the reverse course, which 
is obviously much more difficult to implement. Unfortunately, and in contrast with the 
power spectrum of the primary anisotropy, the shape of the S-Z power spectrum is rather 
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Fig. 5. - Cumulative number counts in redshift space as a function of as- The curves correspond 
to the same us values as in fig. 4. 



featureless, implying a high degree of model and parameter degeneracy. This means that 
there are probably several parameter and model combinations that result in the same 
shape for the power spectrum. Performing a best likelihood test would necessitate a large 
collection of parameters characterizing both the background cosmology (e.g. ^lm^ h), 
the large scale structure (cts), and cluster properties {e.g.lC gas morphology, tempera- 
ture, evolution) . The situation will improve significantly when high quality, resolved S-Z 
and X-ray measurements of individual clusters are available, along with more precisely 
determined cosmological parameters deduced from independent CMB experiments and 
large scale structure surveys. Various aspects of the task of removing parameter degen- 
eracies in the analysis of S-Z and X-ray databases have been explored. Mei & Bartlett 
(2003) have studied the possibility of removing the notorious Q.m - erg degeneracy by com- 
bining the angular correlation function of S-Z clusters with the observed local abundance 
of X-ray clusters. Diego & Majumdar (2004) have shown how the observed S-Z power 
spectrum and number counts can be combined to construct a "hybrid power spectrum" 
whose advantage over the conventional power spectrum is its milder sensitivity to (Ts and 
the basic cosmological parameters, thereby increasing its diagnostic power of intrinsic 
cluster properties. Needless to say, there is still much room for further theoretical study 
of S-Z power spectrum and number counts in order to devise optimal methods of analysis 
of the wealth of high quality S-Z measurements expected in the near future. 

4. — Measurements 

In the first two decades following the discovery of the S-Z effect measurements were 
made with single dish ground-based telescopes. These observations (which were reviewed 
in Birkinshaw 1999) are plagued by the need to account for the fluctuating atmospheric 
emission. Observational and modeling uncertainties are particularly bothersome when 
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Fig. 6. - S-Z image of the cluster MS1054 obtained by observations at 28.5 GHz with the 
interferometric BIMA array (from Carlstrom et al.2002). Colors mark the magnitude of the 
(negative) intensity change, which is highest (red) in the cluster center. 



observations are made at a single frequency. Multi-frequency measurements are essen- 
tial for separating out the various contributing signals, including atmospheric emission, 
emission from Galactic dust, cluster radio sources, and CMB anisotropy. 

Most observational work during the last decade has been made with interferometric 
arrays, whose major advantages over single dish telescopes are increased sensitivity to 
specific angular scales and to signals which are correlated between array elements, insen- 
sitivity to changes in atmospheric emission, and high angular resolution on small scales 
(which is needed to subtract signals from discrete radio sources). The improved sensi- 
tivity of radio receivers made it feasible - mainly through the use of low-noise HEMT 
amplifiers - to image the effect in moderately distant clusters, first with the Ryle tele- 
scope (Jones et al. 1993), and then mostly with the BIMA and OVRO arrays (Carlstrom 
et al. 1996, 2001). Extensive observations (reviewed by Carlstrom et al.2002) with the 
latter two arrays, operating at frequencies ^ 30 GHz, yielded high S/N images of some 
60 clusters in the redshift range 0.17 < z < 0.89, one of which is shown in fig. 6. 

Interferometric S-Z observations were made also with the CBI, an array of 13 small 
(0.9 m) dishes, with 3' — 10' resolution, operating in the 26 — 36 GHz spectral range. In 
contrast to the BIMA and OVRO arrays, the compact configuration of the CBI makes 
this imager suitable for observations of nearby clusters. Results from measurements of 
the effect in 7 nearby clusters were reported by Udomprasert et al. (2004). 

Multi-frequency measurements of the effect (on the R- J and Wien sides) are, of course, 
highly desirable for both the extraction of additional information from its spectral shape, 
and for separating the signal out from the various sources of confusion. The SuZIE 
experiment is the first multi-element array operating at high frequencies. Some fifteen 
clusters were observed at three frequency bands (centered on 150, 220, and 350 GHz) 
with the (two versions of the) experiment (Holzapfel et al. 1997a, 1997b, Mauskopf et 
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al. 2000, Benson et al. 2004). High frequency measurements were also made with the 
balloon-borne PRONAOS (Lamarre et al. 1998) - the first detection of the S-Z effect 
from the stratosphere - and groimd-bascd MITO (Dc Pctris ct al. 2002) experiments. 
Finally, high spatial resolution measurements of the effect in the very luminous cluster 
RX J13471145 were made with the NOBA and SCUBA bolometer arrays (Kitayama et 
al. 2004). 

5. — Results from S-Z & X-ray Measurements 

The main virtues of the S-Z effect that make it a uniquely important cosmological 
probe are its physically well understood origin and its (essentially) redshift independence. 
High sensitivity measurements of individual clusters directly yield the integrated pressure 
of hot IC gas, and thereby also the total mass of the cluster. In general, IC gas density 
and temperature profiles outside the cluster central region can be more directly obtained 
from spatially resolved S-Z measurements of A/((r) than from measurements of the (more 
steeply falling) X-ray surface brightness profile. When S-Z measurements are made close 
to the crossover frequency, the cluster velocity along the los can also be deduced. 

The angular diameter distance, and therefore Hq, can be determined from S-Z and 
X-ray measurements. This SZX method to measure Hq has clear advantages over the 
traditional galactic distance ladder method; these include its well understood physical 
nature, and the capability to test that the global expansion is indeed isotropic. It is also 
possible to determine the principal contributions to - flm, and the currently favored, 
Qa - from the Hubble diagram and from the mean value of the gas mass fraction measured 
in a (sufficiently large) sample of clusters. 

The feasibility of detecting clusters at large redshifts strongly motivates carrying out 
number counts through cluster surveys in order to characterize the population and its 
cosmological evolution. The power spectrum of the CMB anisotropy induced by clusters 
can yield important information on the cluster mass function, cluster properties, and 
the evolution of clusters. Finally, the anisotropy and redshift evolution of the CMB 
temperature, T(z), can be determined from multi-frequency measurements of the effect 
in clusters at different sky directions and redshifts. 

In the following subsections we review some of the recent results obtained from current 
S-Z measurements. The discussion here is very brief and devoid of any considerations of 
important observational issues, such as impact of confusing signals and related systematic 
uncertainties. Observational aspects are expertly reviewed and assessed by Birkinshaw 
in this volume. Comprehensive discussions of results from observational S-Z work can be 
foimd in the reviews by Birkinshaw (1999) and Carlstrom et al. (2002). 

The gas mass fraction fgi The los integrated gas density can be directly determined 
if the projected temperature profile is known from spatially resolved spectral X-ray mea- 
surements. For a relaxed cluster the measiircd gas density and temperature can be used 
in the hydrostatic equation to deduce the total cluster mass M(r) interior to a radial po- 
sition r. The ratio of the gas to total cluster mass, fg, is a good measure of the baryonic 
mass fraction of the cluster, /{,. Determining this ratio for a sufficiently large region, such 
as rsoo - the radius at which the mean mass density of the cluster is 500 times the back- 
ground value - (not only well samples this fraction in the cluster but also) provides a fair 
estimate of the universal value of the baryonic density parameter (e.g. Evrard 1997) when 
averaged over a sufficiently large sample of clusters. Current best S-Z estimate for fg is 
based on BIMA and OVRO measurements of 18 intermediate distance (0.14 < z < 0.83) 
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clusters, which yield fgh ~ (0.08 ± 0.01)/i~^ for the currently favored matter, = 0.3, 
and cosmological constant, CIa = 0.7, density parameters (Grego et al. 2001). The de- 
pendence of the gas fraction on the cluster angular diameter distance can be used to 
determine flm- Analysis of the above dataset yielded the estimate ~ 0.25 in a flat A 
dominated model with h = 0.7. 

Hq and O: The Hubble constant can be determined from measurement of the angular 

diameter distance of a cluster, cIa, from resolved S-Z and X-ray observations. When 
the observed clusters span a wide redshift range, the main contributions to the density 
parameter - Sim and Qa - can then be deduced from the Hubble diagram (a plot of ^^(z)). 
Analysis of a database of 38 measured distances to 26 clusters (from single dish and 
interferometric BIMA & OVRO measurements) at redshifts z < 0.83 yielded a mean value 
of Ho = 60 ± 3 kms~^ Mpc~^ for a flat model with Om = 0.3 and Oa = 0.7 (Carlstrom 
et al. 2002). This ^ 5% {la) observational error is much smaller than the estimated 
systematic uncertainty of ~ 30%. Contributions to the latter include the unknown gas 
thermal profile, the assumptions of sphericity of the gas spatial configuration, possible 
small scale clumping, confusion due to CMB primary anisotropy, and the (unknown) 
cluster peculiar motion (see Birkinshaw 1999 for details). 

Clearly, in order for this method for the measurement of Hq to yield a value whose 
level of precision is comparable to those obtained by the local distance ladder, the CMB 
anisotropy, and SN la methods, the overall level of uncertainty has to be lowered to ~ 5%. 
This could be attained when high quality resolved measurements of a large sample of 
clusters are made. In order to achieve optimal level of overall systematic uncertainty, the 
sample must include many nearby clusters. 

Cluster velocities: Measurements of cluster radial velocities (in the CMB frame) from 
the kinematic S-Z component necessitate sensitive observations in a narrow spectral band 
near the crossover frequency (where the thermal component vanishes) . The exact spectral 
shape of the thermal component near this frequency needs to be known. Confusing signals 
due to the CMB primary anisotropy may very well be the main limitation in reaching 
levels of precision of a few hundred km/s. Attempts to measure radial velocities using this 
method were made mostly with the SuZIE experiment (Holzapfel et al. 1997b, Mauskopf 
et al. 2000, Benson et al. 2003). The SuZIE radial velocity sample includes 8 clusters with 
velocities in a wide range, with uncertainties so large that none of the deduced values 
are statistically significant. 

CMB temperature: Measurements with the COBE/FIRAS experiment have shown 
that the CMB spectrum is a precise Planckian with Tq = 2.725 ± 0.002 K at the current 

epoch (Mather et al. 1999). In the standard cosmological model, T{z) = To(l + z), a 
fundamental relation which has not yet been fully confirmed observationally. Cosmologi- 
cal models with a purely blackbody spectrum but with a different T{z) dependence than 
in the standard model are - formally, at least - unconstrained by the FIRAS measure- 
ments. Also unconstrained are models with spectral distortions that are now negligible, 
but may have been appreciable in the past. Thus far T{z) has been determined mainly 
from measurements of microwave transitions in interstellar clouds in which atoms and 
molecules are excited by the CMB {e.g. LoSecco et al. 2001). The temperature was deter- 
mined in the Galaxy, as well as in clouds at redshifts up to 2: ~ 3. Results are, however, 
substantially uncertain due to the poorly known physical conditions in the absorbing 
clouds. 

The use of the thermal S-Z effect to measure T{z) was suggested long ago (Fabbri, 
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Melchiorri & Natale 1978, Rephaeli 1980). The method proposed by Rephaeli is based 
on the steep frequency dependence of A/j on the Wien side, and the weak dependence 
of ratios of the intensity change at different frequencies on the properties of the clus- 
ter. FormaUy, in the nonrelativistic limit such a ratio is completely independent of the 
Comptonization parameter. Most of the dependence on the cluster parameters drops 
out also in the exact rclativistic description, but a weak dependence remains on the gas 
temperature; the (unknown) cluster velocity introduces a small systematic uncertainty. 
S-Z measurements have the potential of yielding more precise values of T{z) than can be 
obtained from ratios of atomic and molecular lines. 

The availability of spectral measurements of the S-Z effect enabled implementation of 
the method of Rephaeli (1980) to measure T{z) in the Coma and A2163 clusters; the mea- 
surements and their analysis arc described by BattistcUi ct al. (2002), and by Melchiorri 
& Olivo- Melchiorri in this volume. Spectral measurements of Coma {z = 0.0231 ±0.0017) 
and A2163 {z = 0.203 ± 0.002) at four frequency bands yield three independent intensity 
ratios for each cluster; all combinations of these ratios were compared to the theoretically 
predicted values. Fits of the measured ratios to the predicted values were performed, 
yielding best fit values for the CMB temperature at the redshifts of the two clusters, 
Tcorna = 2.789^0 065 ^ and Ta2163 = 3.377^0 103 K (at 68% confidence). These values 
are consistent with those expected from the standard relation T(z) = To(l + z). Battis- 
telli et al. (2002) have also tested two alternative scaling relations that are conjectured 
in non-standard cosmologies, T{z) = To{l + z)^"", and T{z) = To[l -|- (1 -|- d)z] (e.g., 
Lima et al. 2000). They determined the best fit values for the two parameters to be 
a = — 0.16to!32i ^^'i = 0.17 ± 0.36 (at 95% confidence), values that are consistent 
with zero, so no significant deviation was found from the standard model. LoSecco et 
al. (2001) obtained a = -0.05±0.13 and d = 0.10±0.28 (at 95% CL) from measurements 
of microwave transitions. The two sets of results arc consistent. Thus, the S-Z results of 
Battistelli et al. (2002) already provide the same level of precision even though the two 
clusters are at much lower redshifts than the galaxies in the sample used by LoSecco et 
al. (2001). With more precise spectral S-Z measurements expected in the future, it is 
anticipated that the S-Z method will provide a preferred alternative to the atomic and 
molecular lines method. 

6. — Prospects for the Near Future 

The quality of the scientific yield from the many S-Z images obtained with the in- 
terferometric BIMA and OVRO arrays proved beyond doubt that the S-Z effect is an 
indispensable cosmological probe whose great potential has just begun to be exploited. 
New projects have either already begun observations of the S-Z effect or will become 
operational in the next few years. These include the ground based AMI, ACT, AMiBA, 
APEX, MAD/MITO, and SPT, the stratospheric project OLIMPO, and the Planck satel- 
lite, whose full sky survey is expected to result in the detection of thousands of clusters 
(as reviewed by Hansen in this volume). 

With the much improved multi-frequency and high spatial resolution capabilities of 
essentially all the new S-Z projects, the effect will be mapped in many hundreds of 
clusters, in both pointed and survey modes, thereby greatly enhancing the scope of the 
measurements and the quality of the scientific yield. Deep, high resolution measurements 
of the effect in many nearby clusters will provide the best database for precise charac- 
terization of the S-Z properties of clusters. The detailed spectral and spatial images will 
not only yield state-of-the-art gas and total mass profiles, but will also provide us with 
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the key knowledge to assess and reduce the overall level of systematic uncertainties. This 
will lead to a great improvement in the precision of the derived values of cluster masses 
and of the Hubble constant. Extensive cluster surveys and the measurement of the S-Z 
induced CMB anisotropy will yield the mass and redshift distributions of clusters. Im- 
portant information will be extracted on the evolution of clusters and the large scale 
structure. Cosmological parameters, including the dark energy content of the universe, 
will be determined at a high level of precision. 
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